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We review recent advances in imaging the flow of concentrated suspensions, focussing on the 
use of confocal microscopy to obtain time-resolved information on the single-particle level in these 
systems. After motivating the need for quantitative (confocal) imaging in suspension rheology, 
we briefly describe the particles, sample environments, microscopy tools and analysis algorithms 
needed to perform this kind of experiments. The second part of the review focusses on microscopic 
aspects of the flow of concentrated model hard-sphere- like suspensions, and the relation to non-linear 
rheological phenomena such as yielding, shear localization, wall slip and shear-induced ordering. 
Both Brownian and non-Brownian systems will be described. We show how quantitative imaging 
can improve our understanding of the connection between microscopic dynamics and bulk flow. 
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I. INTRODUCTION 

Understanding the deformation and flow, or rheology, 
of complex fluids in terms of their constituents (colloids, 
polymers, or surfactants) poses deep fundamental chal- 
lenges, and has wide applications [1]. Compared to the 
level of understanding now available for the rheology of 
polymer melts, the rheology of concentrated suspensions 
lags considerably behind. In the case of polymer melts, 
it is now possible to predict with some confidence flow 
fields in complex geometries starting from a knowledge of 
molecular properties 0. No corresponding fundamental 
understanding is yet available for concentrated colloids. 
The reasons for this are as follows , Figure [TJ 

In a polymer melt, each chain moves under the topo- 
logical constraints imposed by many other chains. The 
number of these constraints, typically of order 10 3 , is 
sufficiently large that a mean field picture, the so-called 
tube model, can be successfully applied [l(. Moreover, 
the topological entanglement between chains means that 
the breaking of covalent bonds is needed to impose large 
deformations, so that strains often remain small (<C 1). 
In contrast, the maximum number of neighbours in a 
(monodisperse) suspension of spheres is about 10, so that 
mean-field averaging of nearest-neighbour cages will not 
work, and local processes showing large spatio-temporal 
heterogeneities are expected to be important. More- 
over, no topological constraints prevent the occurrence 
of strains of order unity or higher, so that very large 
deformations are routinely encountered. These two char- 
acteristics alone render suspension rheology much more 
difficult. Added on top of these difficulties is the fact that 
concentrated suspensions are in general non-ergodic, i.e. 
they are stuck in some solid-like, non-equilibrium amor- 
phous state. On the other hand, in a monodisperse sys- 
tem or in a mixture with carefully chosen size ratios, 
highly-ordered (crystalline) states can occur. In cither 
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FIG. 1: A schematic comparison between polymer melt rhe- 
ology and colloid rheology. (a) In a polymer melt, a typi- 
cal chain (dashed curve) is constrained by many (in reality, 
~ 10 ) other chains, here represented by small circles. This 
gives rise to the fruitful mean-field concept of a tube in which 
the chain has to move, (b) In a concentrated colloidal sus- 
pension, a typical particle (hatched) is surrounded (in 3D) 
by ~ 10 neighbours. This number is too small for mean-field 
averaging to be meaningful, (c) Large deformations in poly- 
mer melts, such as the process (i) — > (ii), involves breaking 
covalent bonds, and so do not ordinarily occur, (d) There are 
no covalent constraints on order unity deformations, such as 
(i) — > (ii), in a colloidal suspension. 



case, any flow necessarily entails non-linearity (yielding, 
etc.). Moreover, a suspension is a multiphase system, so 
that the relative flow of particles and solvent can, and 
often does, become important. This complication does 
not arise in polymer melts. Finally, compared to the 
considerable effort devoted to the synthesis and experi- 
mental study of very well characterized model materials 
in polymer melts, corresponding work in colloid rheology 
remains relatively rare. It is therefore not surprising that 
the rheology of concentrated suspensions is not nearly as 
well understood as that of polymer melts. 

Advances in this field will pay rich dividends - the suc- 
cessful processing and application of concentrated col- 
loids more often than not depends on understanding, tun- 
ing and exploiting their unique flow properties But 
to build and validate predictive theories of bulk rheolog- 
ical properties, we need microstructural information. A 
case study of this claim comes from recent work on the 
qualitatively distinct rheologies of so-called repulsive and 
attractive colloidal glasses. 

Colloidal glasses are concentrated suspensions in which 
long-range diffusion effectively vanishes @. These 
dynamically-arrested amorphous states have finite shear 
moduli in the low-frequency limit. In slightly polydis- 
perse hard-sphere like suspensions, the transition from 
an ergodic fluid (where long-range diffusion is possible) 
to aglass occurs at a particle volume fraction of <j> « 0.58 
0j- The cause of dynamical arrest is crowding. As (f> 
increases, each particle spends longer and longer being 
caged by its nearest neighbours, until at <f> « 0.58, the 
life time of these cages becomes longer than any reason- 
able experimental time window. Each particle can still 
undergo Brownian motion within its nearest-neighbour 
cage, but its root mean-square displacement saturates at 




FIG. 2: Cartoon of two different yielding mechanisms in an 
attractive glass under oscillatory shear (motion indicated by 
arrows), (a) Un-sheared configuration of particles at time 
zero. The shaded particles bonding neighbours are particles 
1, 2 and 6 and topological neighbours are particles 1-6. At 
small shear amplitude 70 the particles retain the same topol- 
ogy during shear, (b) Above the lower yield strain, bonds 
break, (c) After one period of oscillatory shear, the shaded 
particle retains the same topological neighbours but differ- 
ent bonding neighbours, cf. (a), (d) If the strain amplitude 
exceeds the higher yield strain, the identity of the topolog- 
ical neighbours around the shaded particle is changed after 
one period of shear, cf. (a), the glass is melted and becomes 
liquid-like. 

just over 0.1a (where a is the particle radius) Q, this 
quantity being a measure of the time-averaged fluctuat- 
ing cage size. When oscillatory shear strain is applied to 
a hard-sphere colloidal glass, yielding occurs in a single- 
step process at a strain amplitude of just over 10% 
This has been interpreted in terms of cage breaking: the 
glass yields when nearest-neighbour cages are strained 
beyond their natural (thermally induced) deformation. 

When a strong enough short-range attraction is present 
(where short means a few per cent of a, a second type of 
glassy state occurs - the attractive glass, in which arrest 
is due to particles being trapped by nearest neighbour 
bonds [9J. It turns out that under oscillatory strain, an 
attractive glass yields in a two-stepped process [To[ • The 
first step occurs at a strain amplitude of a few per cent, 
matching the range of the interparticle attraction, and 
then at a second step when the strain amplitude reaches 
a few tens of per cent. The first step has been inter- 
preted as the breaking of interparticle bonds, Figure [5] 
Once these are broken many times, however, each par- 
ticle realises that it is still in a topological cage formed 
by its neighbours. A second, cage breaking, step is still 
necessary for complete yielding to occur. This two step 
signature is discernible in other rheological tests applied 
to attractive colloidal glasses [11| . 

For the purposes of this review, the important point 
is that the explanation proposed for the observed quali- 
tative difference between repulsive and attractive glass 
rheology have been couched in microstructural terms, 
Figure [5] Thus, conclusive validation of these explana- 
tions depends on studying yielding at the single-particle 
level. Such investigation can, in principle, be performed 
using computer simulations. But simulating particles 
under flow with realistic system sizes and including hy- 
drodynamic interactions is a major challenge, although 
progress is being made (see the useful introduction in [12j 
and references therein). Experimentally, the situation is 
also difficult. Until comparatively recently, the method 
of choice for obtaining microstructural information in col- 
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loids is scattering (X ray, neutron or light), whether in 
quiescent or flowing systems. But scattering can provides 
information on the average structure; obtaining local in- 
formation is difficult, if not impossible. The ideal method 
for gaining local information is direct imaging. 

The appropriate imaging method of imaging colloids 
depends, of course, on a number of factors, principally 
perhaps the size and concentration of the particles con- 
cerned. The size imposes a basic constraint, so that, for 
example, visible light can only be used to image parti- 
cles (at least in the 'far field') with radius > 300 nm. 
The concentration is also an important factor, because 
any multiple scattering by concentrated systems will de- 
grade image quality. In this review, we focus on the use 
of confocal optical microscopy to study colloidal systems 
under flow. Confocal imaging rejects out-of-focus infor- 
mation using suitably placed pinholes, thus permitting 
the study of concentrated systems. While the use of con- 
focal microscopy to study quiescent concentrated colloids 
is, by now, well established [H, HI, [TH, [l6[ , the extension 
to flow is relatively recent. The quantitative use of this 
methodology depends on developments in both hardware 
and software. First, to minimise image distortion, fast 
confocal scanning is necessary; we review recent devel- 
opments in this area in Section [II B II Moreover, sample 
geometries are required for simultaneous flow/rheometry 
and imaging, which we review in Section III CI In terms 
of software, algorithms must be developed to track par- 
ticles from image sequences distorted by flow; such al- 
gorithms are explained in Section IIIII After discussing 
these methodological issues, we turn to review a number 
of applications (Section IIVI) . 

To date, the use of quantitative confocal microscopy to 
image concentrated colloids under flow has been applied 
mostly to the study of very well characterised, model 
systems, particularly systems of particles that interact 
as more or less perfect hard spheres. The focus on model 
systems is partly dictated by the desire to obtain as high 
quality images as possible, which requires the matching of 
the refractive indices of the particles and the suspending 
medium (see Section III A|) . The degree of control neces- 
sary, both in the particle synthesis and in solvent choice, 
is only really achievable in model systems. The reason 
for choosing to study hard spheres is that these repre- 
sent the simplest possible (classical) strongly-interacting 
particles. The properties of a quiescent collection of 
hard spheres are well understood in statistical mechanical 
terms, much of this understanding having been obtained 
in the last few decades through the study of model hard- 
sphere colloids. It is therefore natural to extend such 
studies into the area of driven systems. 

We first briefly summarize the quiescent behaviour of 
ideal hard spheres. As first shown via numerical simu- 
lations [It], EH, below a volume fraction of 4>f — 0.494, 
the thermodynamically stable phase is a fluid of colloidal 
particles. Long-ranged order is absent and the particles 
are able to explore all available space (the system is er- 
godic). For 0.494 < (j> < 0.545, the equilibrium state 



is a coexistence of a fluid phase and a crystalline phase. 
At cf>c — 0.545 the entire system must be crystalline to 
minimise the free energy, and this remains the case up to 
the (crystalline) closed packing fraction of 0cp = 0.74. 
This behaviour was largely confirmed by the experiments 
of Pusey and van Megen Q using sterically-stabilised 
PMMA particles (for which see the next section). But, 
they also found that their system failed to crystallize for 
volume fractions <j> > 4>q = 0.58, remaining 'stuck' in a 
non-ergodic (or glassy) state up to the maximum pos- 
sible concentration for amorphous packing (which, for 
monodisperse hard spheres, is </>rcp — 0.64). We have 
already introduced such colloidal glasses. Our review of 
microscopic phenomena will cover both the flow of amor- 
phous states (Section IIV A[) and of ordered (or ordering) 
states (Section |IVB[) of concentrated colloids. While our 
focus will be on hard-sphere colloids, we will also mention 
various imaging studies on the the flow of non-Brownian 
suspensions. 



II. EXPERIMENTAL METHODS 
A. The colloidal particles 

One of the key ingredients for high-quality confocal 
imaging of concentrated suspensions is a colloidal sys- 
tem that allows identification of the individual spheres, 
does not suffer from rapid photobleaching and gives 
good depth penetration into the bulk of the disper- 
sion. Fluorescently-labelled poly-methyl-methacrylate 
(PMMA) particles in hydrocarbon solvents show these 
characteristics. 

The synthesis of har d-sp here PMMA particles was first 
described by Barrett [lj| in 1975 and subsequently by 
Antl et al. [2(| in 1986. It is a two-stepped disper- 
sion polymerization reaction, yielding particles made of 
PMMA cores kept stable by a thin (~10nm) outer layer of 
poly-12-hydroxystearic acid chains, which act as a steric 
barrier to aggregation. In the first step of the polymer- 
ization the spheres are made by growing PMMA chains 
in solution, which become insoluble when they reach a 
certain size. At this point they come out of solution 
and clump together to form particles that are kept sta- 
ble by physisorbed poly-12-hydroxystearic acid chains. 
The second step of the preparation involves chemically 
linking the stabiliser chains to the spheres. These chains 
are tightly packed on the particle surface [2l[ and stretch 
out in good solvents such as various hydrocarbons, which 
causes the particles to interact as nearly-perfect hard- 
spheres [22j]. The resulting particles can have polydis- 
persities as low as a few percent, and can be made with 
radii ranging from 100 nm up to and above 1 /im. 

Fluorescent labelling of the particles for confocal mi- 
croscopy may be achieved in three ways. The first in- 
volves the use of polymerisable dyes. These dyes have 
been chemically modified to include a reactive group 
that can be chemically attached to the particle as they 
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are produced. The advantage of this procedure is that 
the dye will not leave the particle once it is incorpo- 
rated. For sterically-stabilised PMMA particles this in- 
volves adding a methacrylate group to the dye, and 
several such procedures [H, [H, [2j| [H, [2?| have been 
described in the literature. The most commonly used 
dye is 7-nitrobenzo-2-oxa-l,3-diazole-methyl methacry- 
late (NBD-MMA) [H, H3, which is excited at 488nm 
and emits at 525nm, while the red end of the spectrum is 
well served by (rhodamine isothiocyanate)-aminostyrene 
(RAS) [22|. " 

The second way to dye the PMMA spheres is to add 
an unreactive species during particle formation. Here the 
dye has no polymerisable group and is just dissolved in 
one of the reaction reagents with the hope that it will 
become incorporated into the growing particle. The ad- 
vantage of this method is that no chemistry, which may 
alter the dye's physical properties, is required prior to 
use; it also allows for a wider range of dyes. The dis- 
advantage is that the dye may leak out when, for exam- 
ple, solvency conditions are changed, or migrate within 
the particle. This technique was employed by Camp- 
bell and Bartlett [28[, who examined how four differ- 
ent red dyes affected particle formation. They found 
optimum properties when using l,l-dioctadecyl-3, 3,3,3- 
tetramethylindocarbocyanine (DiIC18), which had sig- 
nificantly slower photobleaching rate than other dyes 
tested, did not affect particle preparation and did not 
interfere with the hard-sphere behaviour. However, this 
dye degrades at temperatures around 100 °C and there- 
fore the reaction that chemically links the stabiliser to 
the spheres can not be performed (29j . 

The third method to stain the PMMA is to add 
the fluorescent dye after particle synthesis. This is 
achieved by finding a solvent that will dissolve the 
dye and also be taken up by the particles. Thus, 
an acetone/cyclohexanone mixture can be used to de- 
liver rhodamine perchlorate dye to preformed PMMA 
spheres [l3[ • The advantage of this method is that once 
a suitable delivery system is found, many possible dyes, 
and even multiple dyes, may be added to the spheres. 
The disadvantage is that the solvent mixture may attack 
the spheres, swell them or alter their physical properties. 

An advantage of PMMA spheres is that the polymer 
may be cross- linked 30] . This is achieved via a molecule 
with two polymerisable groups which is used to chemi- 
cally bind all the individual PMMA polymer chains in 
a particle together into a network. This can stabilise 
the particles in solvents which would normally dissolve 
them such as aromatics. An additional benefit is that 
this method allows the fluorescent dye to be kept within 
one particular area of the particle, usually its core. The 
preparation of such core-shell particles [23, bLfl, [33, [33| in- 
volves modifying the usual reaction by adding the cross- 
linking agent and dye at the start of the procedure and 
creating particles as usual. However, instead of going 
on to chemically link the stabiliser to the particles, they 
are cleaned to remove excess dye and then more methyl 



methacrylate monomer and cross-linking agent are re- 
acted with them to produce an undyed shell. The sta- 
biliser is then chemically attached. Such core-shell par- 
ticles allow for more accurate detection of the particle 
centers from microscopy in concentrated systems, as the 
fluorescent cores are well separated from each other (e.g. 
0). 

To achieve good imaging conditions in the bulk of the 
sample and reduce scattering of the laser and excited 
light, the refractive index (RI) of the particles and the 
solvent should be closely matched. For PMMA spheres 
the RI is around 1.5 and can be matched using a mixture 
of cis-decahydronaphthalcnc (RI=1.48) and tetrahydron- 
aphthalcnc (RI=1.54) in a ratio of approximately 2:1. In 
earlier work [6j a decalin-carbon disulphidc mixture was 
used (ratio 2.66:1) but this is problematic due to the 
volatility and toxicity of the carbon disulphide. The par- 
ticle RI may also be modified by a few per cent by adding 
different monomers during preparation [35| . 

Another concern when studying colloidal dynam- 
ics is sedimentation. To counter this phenomenon, 
solvent mixtures have been sought which not only 
index-match the colloids but also allow density 
matching. Adding carbon tetrachloride (36| to the 
cis-decahydronaphthalene /tetrahydronaphthalcnc sys- 
tem was tried but had only limited success as it 
enhanced photobleaching and imparted a charge to the 
particles. Another mixture used consists of cis-decalin 
and cycloheptylbromide (CHB) (3|| (3?J (or cyclohexyl- 
bromide (38|). This also imparts a charge to the colloids, 
likely due to photo-induced cleavage of the Br-C bond 

and subsequent solvent acidification, but this may 

be screened by the addition of salt [3^, |4(j, (partly) 
restoring the nearly-hard-sphere nature of the system. 

For both RI and density matching, it is important to 
note that the extra solvent components can swell the 
PMMA particles, and the solvent uptake may take up 
to weeks to saturate. Time-dependent monitoring is the 
only sure means of ensuring that the particle size has sta- 
bilised. Otherwise, significant errors in volume fraction 
estimation may result. 

Besides PMMA particles, various other systems have 
been considered in the literature for use in confocal mi- 
croscopy. One of the best known is silica spheres. They 
are prepared by the hydrolysis and condensation of alkyl 
silicates in ethanol using ammonia as catalyst as first 
described by Stober et al. [4l[. They can be made flu- 
orescent by adding dyes that have been reacted with 
silane coupling agents which make them affix to the sil- 
ica [421 . 143 . [4J, [4^| . The advantage of this system is 
that various dyes can be used and added at any time 
during the synthesis [12]. The spheres are stabilized in 
organic solvents via a dense la yer of organophilic material 
grafted onto their surface [H, H3, IH, [4t| , and the RI can 
be matched with the solvent (RI=1.45 for Stober silica 
spheres). However, the main problem with silica spheres 
is that they are quite dense (rep orted values range from 
1.51g/cm 3 49] to 2.2g/cm 3 [501] ) so that sedimentation 
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problems may be severe with all but the smallest spheres. 

Water-based fluorescent particles can also be used, but 
their RI differs considerably from that of water, limiting 
observation in the bulk. The RI of water can be modi- 
fied by adding salts [5l| , but the concentrations required 
are exceptionally high, making most particles unstable 
against van der Waals attraction. 

A final interesting development is that of quantum-dot- 
loaded particles. Quantum dots [5j| are themselves very 
small semiconducting colloids (1 nm to 10 nm) which 
fluoresce due to quantum confinement. They can be 
trapped within a polymer or silica [531 ] colloid. Since 
quantum dots photobleach much less than organic dyes, 
they can be used for experiments involving long-term ob- 
servation. 



B. Imaging 

Direct imaging dates back to the work of Robert Brown 
|54| , who used it to discover and study Brownian motion, 
which is the defining characteristic of colloids. Subse- 
quently, Perrin [55| used direct imaging to great effect 
in his Noble-Prize-winning work on sedimentation equi- 
librium and diffusion in dilute suspensions. However, di- 
rect imaging has flourished as a major research tool in 
colloidal research only in the last two or three decades, 
mainly due to the increase in imaging and computing 
power. 

Our main emphasis in this article is on high-resolution 
imaging and reconstruction of the position of every par- 
ticle in some volume of a concentrated colloidal sus- 
pension under flow in real-time. Before reviewing this 
subject in detail, we point out the use of conventional, 
low-resolution imaging to track the position of a num- 
ber of tracer particles. Perhaps most important for our 
purposes here, the well-established technique of parti- 
cle imaging velocimetry (PIV) [56| can be used to shed 
light on complicating factors in conventional rheological 
measurements such as wall slip I57J . [5 81 and flow non- 
uniformities such as shear banding [59j |. Thus this tech- 
nique has recently been used in a rheometer to elucidate 
the physics of wall slip in concentrated emulsions under 
shear [6(| The correlation techniques used in PIV 
to measure tracer velocity can also be used to track the 
shear-induced diffusion of tracers in concentrated suspen- 
sions of non-Brownian spheres [62|, [63| . PIV and related 
techniques are, of course, limited to transparent samples. 

For completeness, we mention other methods for ve- 
locimetry that have no requirement for transparency, 
such as heterodyne light-scattering (64[ and ultrasonic 
velocimetry [65j . The latter has been applied to charac- 
terise slip and flow non-linearities in micelles and emul- 
sions (66l . We also mention under this heading Nu- 
clear Magnetic Resonance Imaging (NMRI) [H, Hil [zO, 
[7lj . another velocimetry method independent of trans- 
parency, which can also provide information on local den- 
sity. The technique has spatial resolution down to ~ 



20 fim and has been combined with rheometric set ups to 



relate velocity profiles to macroscopic rheology 72 , |73j , 
to give insight on the occurrence of shear bands 74[ and 
shear thickening [75| in concentrated suspensions. 

While all of these techniques give additional insight 
unavailable from bulk rheology alone, building up a com- 
plete picture of colloidal flow requires dynamical infor- 
mation on the single-particle level. We now turn to mi- 
croscopic methods that give precisely such information. 
We focus on single-particle imaging in 3D, but mention 
that the imaging of a single layer of colloids has been 
used to great effect to study fundamental processes in 
2D (e.g. [7l,[z3,[2l|)- While perhaps somewhat less com- 
plex than 3D imaging, 2D imaging nonetheless presents 
some challenges, e.g. when the imaged objects come into 
very close proximity [791 ]. 

The use of conventional (non-confocal) optical mi- 
croscopy to study concentrated colloidal suspensions in 
3D has been reviewed before [80]. In nearly index- 
matched suspensions, contrast is generated using either 
phase contrast or differential interference contrast (DIC) 
techniques. One advantage of conventional microscopy 
is speed: image frames can easily be acquired at video 
rate. Conversely it has poor 'optical sectioning' due to 
the presence of significant out-of-focus information, so 
that particle coordinates in concentrated systems cannot 
be reconstructed in general, although structural informa- 
tion is still obtainable under special circumstances 811 ] . 

Compared to conventional microscopy, confocal mi- 
croscopy delivers superior 'optical sectioning' by using 
a pinhole in a plane conjugate with the focal (xy) plane. 
It allows a crisp 3D image to be built from a stack of 2D 
images even for somewhat turbid samples, but each 2D 
image is acquired by scanning, which imposes limits on 
acquisition speed. The technique has been described in 
detail before |2]. 

The use of confocal microscopy to study concentrated 
colloidal suspensions was pioneered by van Blaaderen 
and Wiltzius [HI], who showed that the structure of a 
random-close-packed sediment could be reconstructed at 
the single-particle level. Confocal microscopy of colloidal 
suspensions in the absence of flow has been recently re- 
viewed [HI, [T3, [TH, [Hj]. We refer the reader to these 
reviews for details and references. Here, we simply note 
that this methodology gives direct access to local pro- 
cesses, such as crystal nucleation (§4l [ and dynamic het- 
erogeneities in hard-sphere suspensions near the glass 
transition (§4l [371]. 

Our main interest is the use of microscopy to study 
the flow of concentrated suspensions at single-particle 
resolution in 3D. It is possible to use conventional (non- 
confocal) video microscopy for this purpose (§Bl l§q . [§7l . 
[§§[, but the poor optical sectioning hinders complete, 
quantitative image analysis. Conversely, crisp confocal 
images in principle permit the extraction of particle coor- 
dinates, but due to slow scanning and acquisition rates, 
early observations in real time (i.e. during shear) pro- 
duced blurred images that again limited the potential for 
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quantitative analysis . A common solution was to ap- 
ply shear, and then image immediately after the cessation 
of shear, both in 2D 0, and in 3D [H, M, M, HI • 
In the next subsection we review developments in confo- 
cal microscopy that permit faster acquisition and hence 
time-resolved 3D imaging of particulate systems under 
flow. 



1. Confocal microscopy 

Confocal images are built by scanning a laser beam 
across the field of view and collecting the emitted fluo- 
rescent light through a pinhole (Laser Scanning Confo- 
cal Microscopy, LSCM). Traditionally the laser beam is 
scanned across the specimen by two galvo-mirrors which 
gives maximum acquistion rates of the order of 1 Hz, de- 
pending on image size. Technical advances such as the 
use of resonant galvo-mirrors, spinning (Nipkow) discs 
(possibly extended with an array of micro-lenses) and 
acousto-optic deflectors (AODs) have significantly im- 
proved upon such acquisition rates. 

AODs are crystals which act as diffraction gratings. By 
sending a standing sound wave at radio-frequency across 
the crystal the local index of refraction is changed, creat- 
ing a grating which deflects laser beams passing through 
the crystal. By changing the frequency of the sound 
wave, the diffraction angle is changed and therefore the 
field of view can be scanned extremely rapidly. The main 
problem associated with AODs is that the grating de- 
flects light of different wavelengths by different angles 
and therefore obstructs the fluorescent light travelling 
back along the optical path. This is partly resolved by 
combining an AOD with a galvo-mirror and a slit instead 
of a pinhole. The galvo-mirror positions the beam in one 
direction, the AOD scans in the orthogonal direction and 
the fluorescent signal is detected through the slit. The 
slit slightly reduces the rejection of out-of-focus light and 
causes a slight anisotropy in the in-plane point-spread 
function, but the method offers frame rates > 100 Hz for 
512 x 512 pixels images. 

Another technique that improves the scanning speed 
is the use of a spinning (Nipkow) disc 82]. This solution 
operates by illuminating the sample through an array of 
pinholes printed onto a spinning disc and thus achieving 
the scanning of the sample during the disc's revolution. 
Spinning discs can achieve full scanning of the field of 
view by fractions of a revolution and therefore yield frame 
rates of hundreds of Hz. Disadvantages of this method 
are the fixed size of the pinholes, limiting the use of differ- 
ent objectives for optimum operation, and also the strong 
loss of intensity occurring at the disk. Recently, the last 
disadvantage has mostly been overcome by using laser 
illumination in combination with an additional array of 
micro-lenses, strongly increasing theefficiency and con- 
siderably reducing photo-bleaching [95| . 

Regardless of the scanning principle, the confocal mi- 
croscope can be operated either in a 2D mode, i.e. cap- 




FIG. 3: Reprinted from [96J. Projection of a raw 3D image 
stack of 1.7 fim PMMA particles. The image volume, ixt/x 
z ~ 29 x 29 x 15 /j,m 3 (256 x 256 x 76 voxels), was scanned in 
~ 1 s. 



turing time sequences of images at a fixed focal depth z, 
or in a 3D mode, i.e. capturing image stacks obtained 
by scanning the sample volume, Figure [3) For rapid 3D 
acquisition, the best method of scanning along z is to 
use a piezo-element for focus control of the objective. 
The two operation modes impose different limits on the 
acquisition rates for successful quantitative imaging (see 
Section HHB] and Q). 



C. Flow geometries 

Microscope studies of flowing particulate suspensions, 
and soft matter systems in general, require well-defined 
flow geometries to facilitate data interpretation and al- 
low, in principle, for a mapping of the observed dynam- 
ics to the rheological properties of the system. To date, 
many studies have employed flow cells which impose the 
deformation (rate) but do not have the ability to mea- 
sure directly the (shear) stress. Recently this limitation 
has been overcome by combining confocal imaging with a 
rheometer and by inferring (indirectly) stress from pres- 
sure drops in channel flows. 

In all geometries, wall properties play a crucial role in 
the application of shear. The most direct manifestation 
of wall effects is that of (apparent) slip in many complex 
fluids, in particular colloid gels or pastes driven along 
smooth surfaces [13, HH, H3| • While the physics of slip in 
suspensions and pastes is an interesting topic in its own 
right (see Section IIV A 1[) , in most flow geometries the 
practical goal is to minimise slip and transfer shear to 
the suspension. The remedy for wall slip differs between 
various systems; for moderate to very dense particle sus- 
pensions, a coating consisting of a sintered monolayer of 
similar size particles generally provides stick boundary 
conditions. 
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1. Parallel plate shear cells 

Planar shear, or planar Couette flow, is simply im- 
plemented by placing the material between two parallel 
plates much larger than their separation z gap , and trans- 
lating them relative to each other. This can be achieved 
either by fixing one of the plates and moving the other 
one, or by moving them in opposite directions so that, in 
the laboratory frame, the zero velocity plane is situated 
somewhere between the two plates. For a Newtonian 
fluid without slip and sufficiently far form the edges of 
the plates (a few z gap ), the shear rate in the fluid is con- 
stant, 7 = {vt — vb)/ Zg a p, with vt (vb) the top (bottom) 
plate velocity. In the following, we denote the shear ve- 
locity direction by x, the vorticity direction by y and the 
velocity gradient direction by z. 

By construction, parallel plate shear cells can only 
achieve finite strains after which the direction of motion 
must be reversed; they are thus particularly suited for 
oscillatory strain studies. However, in some designs, ex- 
tremely large strains are possible (~ 20 — 50) so that, 
with constant drive velocity, steady shear is effectively 
achieved in each half-cycle. Using a microscope glass 
(cover) slide as the bottom plate allows imaging of the 
suspension during shear. An important requirement is 
that the plates should be strictly parallel; this condi- 
tion is particularly stringent for small plate separations 
e.g. [981 ]. Non-parallel plates induce non-uniform shear 
as well as drifts in the sample due to capillarity effects. 
Other considerations relate to component weight, me- 
chanical resonances, the driving motor and minimising 
sample evaporation. Depending on the desired operation 
range, compromises may be required to optimise either 
long-time stability or high-frequency behaviour. 

Planar shear cells were initially coupled to conventional 
mic r osco pes [H, H|| or light scattering set ups [H, l99l . 
llOOl . Il0l| . With the advance of confocal imaging they are 
now preferentially used in combination with an inverted 
microscope and a confoca l scanner. 

Cohen et al. [l02l 1 1031 ] used a simple design with a 
movable microscope cover slip as lower plate and a fixed 
top plate. The maximum plate separation 100 /tm and 
they were parallel to within 1/tm over the shear zone. The 
lower plate was driven by a piezoelectric actuator with 
displacements up to 90/tm at frequencies < 100 Hz. The 
sample between the plates was in contact with a large 
reservoir of un-sheared bulk suspension. 

A similar design but allowing for larger strains and 
higher shear rates is described by Solomon et al. [9~jj . 
Two tilt goniometers allow the user to tune the paral- 
lelism of the plates; a gap of 150/im is set by a linear 
micrometer. Oscillatory shear was produced by applying 
a sinusoidal displacement with a linear stepper motor. 
The shear rates were in the range 0.01-100 s _1 and strain 
amplitudes in the range 0.05-23. 

A cell optimised for slow shear and large amplitudes 
was developed by the Edinburgh group [96|, Il04l |. The 
operational gap size ranged from ~ 200 — 1000 /mi with 



plates parallel to ±5 /mi over the shear region. The top 
plate was driven at 0.05 — 10 /mis -1 by a mechanical 
actuator with magnetic encoder. The maximum trans- 
lation was L s ~ 1 cm, allowing steady shear up to a 
total strain A7 = L s /Z gap > 1000%. The cell could be 
operated either with the bottom plate fixed or with the 
plates counter-propagating via an adjustable lever sys- 
tem to tune the zero velocity plane. 

A recent, 'state of the art', parallel plate cell is the de- 
sign by Wu et al. [98|, which combines high mechanical 
stability and a modular construction. The plate sepa- 
ration ranges between 20 to 200/tm, and a special pivot 
system was designed to align the plates parallel to the 
highest degree. A piezo-stepper motor provided plate 
velocities ranging from ~ 2 /mis -1 to 10 mms -1 . As 
in the Edinburgh's design, the relative plate motion can 
reach up to ~ 1 cm, allowing for large accumulated strain 
(< 50) and quasi-steady shear. 



2. Rotational geometries 

Application of continuous shear is achieved in rota- 
tional geometries such as cylindrical Couette, plate-plate 
or cone-plate devic es, th e standard geometries used in 
traditional rheology [l05l |. Of these, cone-plate and plate- 
plate geometries are most suitable for microscopic ob- 
servation. Various ways of conventional (microscopic) 
imaging in such geometries have been used, including 
an early direct o bser vation of crystallization in a plate- 
plate rheometer [l06t | and a microscopic study of a con- 
fined, charged suspension under continuous shear in a 
rotational plate-plate set-up [13 ■ Another technique is 
to image tracer particles in the system from the side of 
a cone-plate or plate-plate geometry to obtain the de- 
formation or velo city profile. This has been done both 
qualitatively 107 1 and quantitatively [6lj . 

For completeness we also mention geometries used to 
study non-Brownian suspensions, where the typical di- 
mensions are much larger (particle radius a > 50 /im). 
One is an annular channel formed by two concentric glass 
cylinders wh ere a ring-shaped top plate drives the sus- 
pension jl08| . the other consists of a special truncated 
counter-rotating cone-plate cell [63[. Both set ups allow 
for imaging in multiple directions using different camera 
positions. As an exponent of imagin g wit h multiple cam- 
eras, we mention the recent work in [109j |. where two or- 
thogonally positioned cameras were used simultaneously 
to image tracers in an immersed granular packing sheared 
in a cylindrical Couette cell. Analysis of the data from 
the two cameras provided the full three-dimensional tra- 
jectory of the tracers. 

The first combination of a rotation al shear cell with 
confocal microscopy was described in jllOl llll| . It con- 
sists of a cone-plate, similar to Figure 3Kb), where the 
cone and the plate (a large glass cover slide) could be 
rotated independently. In contrast to Figure [31(b) , the 
objective was located at a fixed radial position where 
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FIG. 4: Reprinted from [9g|. (a) Schematics of the confocal 
rheoscope of the Edinburgh group [§(|. The top arrow marks 
translation of the rheometer head to adjust the geometry gap, 
the horizontal arrow indicates translation of the arm support- 
ing the objective to image at different radial positions r. (b) 
Close up of the central part of the r heos cope, similar to the 
cone-plate imaging system of Derks except that in the 

latter the lower plate can also be rotated, while in the for- 
mer the microscope objective radial position r can be varied, 
(c) Gap profile of a 1° cone-plate geometry, measured in the 
confocal rheoscope with fluorescent particles coated on both 
surfaces. 



the gap height was 1.7 mm. This exceeds the working 
distance of high numerical aperture objectives so that 
only the lower part of the gap could be imaged at single- 
particle resolution, a general limitation of any large-gap 
geometry. The shear rates obtained were in the range of 
10~ 2 to 10 2 s _1 . The height of the zero velocity plane 
could be adjusted by tuning the ratio of the cone and 
plate rotational speed, while keeping the shear rate con- 
stant. Due to the weight of the cone and plate, oscillatory 
shear experiments were limited to low frequencies. The 
main limitation of this set up is the acquisition rate of 
the confocal scanner: particle tracking could be achieved 
only near the zero velocity plane and in systems where 
the out-of-plane motion was slower than the acquisition 
rate, such as s colloidal crystal. 

A more recent development of confocal imaging in a 
'rheometric' geometry, which also permits simultaneous 
measurement of the rheological properties, is the set up 
of Besseling et al. [96[ , see Figure [U Here a rheometer 
was equipped with a custom-built, open base construc- 
tion, providing space for imaging optics. An adjustable 
plate with an imaging slit at the top of the base is covered 
with a circular cover slide forming the bottom surface of 
the geometry. This allows for imaging at various radial 
distances r. The stress-controlled rheometer (AR2000, 
TA instr.) enables all classical rheological tests either in 
plate-plate or cone-plate geometry. In practice a cone 
with an angle of 1° was used, Figure |UJb),(c). The imag- 
ing optics under the plate (piezo-mounted objective to 
scan along z, mirror and lenses) are coupled to a con- 
focal scanner (VTEye, Visitech) to give full 3D imaging 
capabilities as in a standard microscope mount. A major 
advantage is the acquisition rate (200 fps) of the confo- 
cal scanner, which allows 3D imaging of single particle 
dynamics up to rates of ~ 0.1 s^ 1 and measurements 
of in the xy plane up to 100 s _1 . Finally, either the 
cone or the bottom glass plate (or both) can be coated 
with a sintered layer of colloidal particles for the study 
of slip-related problems. 




FIG. 5: A possible sample cell for capillary flow. The capil- 
lary is not drawn to scale. The construction is placed on the 
microscope stage and imaged from below. 



3. Capillaries and micro- channels 

Micro-channels and capillaries are flow geometries that 
occur in many practical applications. They are also inter- 
esting from a fundamental perspe ctiv e, offering insig ht on 
issues such as flow instabilitie s |l!2|. confineme nt jll3l | 
and particle migration effects jll4lll5l Ill6l . Ill7j . When 
the pressure drop over the channel is measured, it is 
possible, in principle, to relate microsco pic o bservations 
to bulk rheological properties (see e.g. [l!8j |). but this 
is non-trivial as it requires steady, uniform flow along 
the c hann el and absence of entrance or confinement ef- 
fects [ng|. 

From the 1970s, a large number of imaging studies of 
non-Brownian suspensions flowing in mm- to cm-sized 
channels have been perform ed via Laser Doppler Ve- 
locimetry flM fT2TI. [l2^. [l23j ] and Nuclear Magnetic Res- 
onance Imaging [HI but very little structural infor- 
mation has been obtained at the single-particle level. Op- 
tical microscopy experiments on channel flows of colloids 
have only recently started to ap pear, often in relation to 
microfluidics applications [l2(jj ]. To avoid image distor- 
tions, channels with square or rectangular cross sections 
are gen erally preferred to cylindrical capillaries. 

Haw [11 21 ] studied the jamming of concentrated hard- 
sphere colloids at the entrance of mm-sized, cylindrical 
capillaries using conventional microscopy. The shape of 
the channels and the method of imaging militated against 
accessing detailed microscopic information. 

Confocal studies in the group of Weeks pllllZi 

used 

rectangular capillaries (50/imx500/im) coupled to a sy- 
ringe as a reservoir to drive suspensions of intermediate 
volume fractions. Due to the high flow rates involved in 
their experiments, particle tracking was impossible and 
instead the flows were studied using image correlation 
techniques and intensity measurements. A very similar 
geometry, but with square channels was used in [l27| to 
study the flow of attractiv e gel s of s ilica particles. 

Isa and co-workers [H, Il28l Il29l | used a slightly dif- 
ferent geometry to study the flow of very dense colloidal 
pastes. Figure [5] shows a sketch of the setup: square or 
rectangular, micron-sized glass capillaries are connected 
to a glass reservoir at one end and flow is driven by suc- 
tion at the other end of the channel. The capillaries 
are either untreated and smooth, or have the inner walls 
coated with a sintered layer of PMMA colloids. 

The recent development of micro-fabrication tech- 
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niques (soft lithography) has opened up the possibility 
of s tudy ing flow in microfluidics geometries. Degre et 
al. jll8l | have performed PIV studies of polymer solu- 
tions seeded with tracers in micro-fabricated geometries 
made by adhering a moulded block of polydimethylsilox- 
ane (PDMS) onto a glass cover slide, obtaining chan- 
nels with a thickness of tens of microns. PDMS has 
many pr actic al advantages, but swells in most organic 
solvents jl3C| . making it incompatible with many col- 
loidal model systems. It also has a relatively low elastic 
modulus causing deformation under high pressures. A 
solution to these problems is to substitute PDMS with 
a photo-curable monomer, as described in [l3l| . It is 
clear that soft lithography offers great flexibility in de- 
signing channels of almost any geometry paving the way 
for high-resolution imaging studies of colloidal flows in 
a variety of complex micro-environments. Such studies 
will be relevant for microfluidic applications, as well as 
for modelling flows in porous and other complex materi- 
als. 




FIG. 6: Reprinted from 96]. Advection profiles in various 
geometries, (a) A uniform 2D shift AX across the entire 
field of view maximises the correlation, (b) 2D case where 
the advected motion is a function of y; the image is then 
decomposed in bins centered at y q , each of which is shifted by 
AX(y q ) to obtain maximum correlation. (C) 3D image where 
the motion is a function of z only. The stack is decomposed 
in slices centered at z r , each of which is shifted by AX(z r ) to 
maximize the correlation, (d) 3D case with y and z dependent 
motion. Decomposition into y and z bins yields the advection 
profile AX(y q ,z r ). 



III. IMAGE ANALYSIS 

To extract the maximum amount of quantitative infor- 
mation from (confocal) images, detailed image analysis is 
necessary. The central components of such analysis are 
particle location and tracking. 



A. General methods 

The basic method to obtain coarse-grained information 
on the flow without details on the dynamics at the par- 
ticle level is Particle Image Velocimctry (PIV) or related 
correlation techniques [56(. Using a sequence of images 
of tracer particles in the suspension or confocal images 
of the full microstructure during flow, a map of the ad- 
vected motion between consecutive 2D images, i — 1 and 
i, or between parts of these images ('tiles'), is obtained 
as the shift, (AX, AY), which maximizes the correlation 
between these images or regions. By repeating the proce- 
dure over a sequence of frames, one obtains the displace- 
ment as function of time, (AX AY(tj)). In general, 
PIV yields a discrete vector field, AR(r Mr ), with dis- 
placement AR in each element r pqr of a 3D image (r pq 
in 2D). In many practical cases however, the flow field 
has a simpler structure, see the examples in Figure [5] 

For 2D images with uniform motion in the xy plane, 
the procedure is applied to the entire image, Figure [SJa), 
giving AR(r,tj) = AX(r,U) = AX(U). For advective 
motion which depends on the position y transverse to the 
flow, Figure [Sfb), the procedure is performed on image 
strips, yielding AX(y q ) discretised at the strip centers 
y q . For 3D images of simple shear, Figure [6fc), the mo- 
tion is a function of z only. The image stacks are then 
decomposed in xy slices at different z and the procedure 
is applied to the individual 2D slices. In channel flow, 



Figure [6}d, the stacks are first decomposed in xy slices 
and then each slice is decomposed in y-bins for which the 
motion is analysed. 

These methods have been used to measure the velocity 
profiles in various dense suspensions. Meeker et al. [60( 
imaged tracers embedded in pastes from the edge of a 
cone-plate geometry. A number of studies have been 
reported in which micro-PIV is applied to the flow in 
channels jl!8l . Il32j | or to obta in displacement profiles in 
parallel plate geometries [l02| ] . 

A related application of the correlati on m ethod has 
been reported by Derks and co-workers [lllj . In their 
study of dense colloids flowing in a cone-plate geome- 
try, a single confocal scan in the velocity-gradient plane 
(y — z plane in their notation, see Figure [7|) was per- 
formed. This was then analysed by shifting and corre- 
lating lines (along y) between consecutive z- values in the 
image, rather than performing PIV between consecutive 
frames. Via this procedure the distortion of the particle 
images could be quantified, yielding parabolic displace- 
ment profiles, Figure [3 corresponding to linear velocity 
profiles from which the shear rate was extracted. 



B. Locating particles 

To go further in the analysis, one first needs to deter- 
mine the location of the particle centers from the images. 
The standard method was developed over a decade ago 
by Crocker and Grier (CG) jl33| . and has since been used 
in numerous studies on colloidal dynamics. 

While some of the concepts in their method ma y als o 
apply to identifying objects with varying shapes [l34| . 
the method is primarily aimed at identifying circular 
(2D) or spherical (3D) objects that appear bright on a 
dark background. Since experimental images have un- 
avoidable pixel noise and often also undesired intensity 
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FIG. 7: Reprinted with permission from 111], copyright 
[2004], Institute of Physics Publishing. Confocal images (yz, 
75 /Ltmx56 /xm, 512x512 pixels 2 ) of a colloidal fluid at vari- 
ous shear cond itions taken in the counter-rotating cone-plate 
shear cell |lll| ] , see Section III C 21 The applied shear rates 
are 1.67, 3.36 and 8.39 s _1 (top to bottom); the ratios of the 
applied cone to plate rotation speeds are 84, 129 and 175 (left 
to right). Graphs A and B show displacement profiles, y(z), 
measured from these images via cross-correlation of scanned 
lines. The appropriate profile is overlaid on each image (white 
curves). The velocity profiles (dy/dz) calculated from these 
displacement profiles are shown in the graphs C and D. The 
particle diameter is 1.50/im. 



modulations, the images are first treated with a spatial 
bandpass filter, which eliminates long-wavelength con- 
trast gradients and pixel-to-pixel noise. Next, the co- 
ordinates of the centers of the features are obtained by 
locating the local intensity maxima in the filtered images. 
These coordinates are then refined to a higher accuracy 
by applying a centroiding algorithm, which locates the 
brightness- weighted center of mass (centroid) of the par- 
ticles; with this refinement the coordinates of the particle 
center can be ob taine d with a typical resolution of 1/n 
of the pixel size [l33j ] where n is the particle diameter 
in pixels. However, the method has some limitations in 
concentrated systems, as individual particle images may 
start to overlap. Unless one uses the core-shell particles 
mentioned in Section III Al an alternative method to the 
coordinate refinement may be required. A very useful 
technique for doing this, based on optimizing the overlap 
between the measured intensity profile of each particle 
and the so called 'sphere spread function' , ha s been de- 
scribed recently by Jenkins and Egelhaaf [l35| |. 

In addition to the intrinsic sub-pixel accuracy from the 
above refinement methods, additional errors arise also 
from particle motion during image acquisition. These er- 
rors may be considerable when images are obtained by 
a scanning system, where the pixels are not acquired in- 
stantaneously. For an acquisition time l// sca n for a 2D 
image with n lines and a particle radius (in pixels) of a, 



the acquisition time of a particle image is: 

C = 2a/(n/ scan ). (1) 

For a 3D image (a z-stack of 2D slices), the voxel size 
in the z direction may differ from that in the x and y 
direction. When the particle radius in z-pixels is given 
by a z , the acquisition time for a 3D particle image is: 

^im = ^^z//scan- (2) 

Typical parameters (/ sca n = 90 Hz in a fast confocal, 
n = 256, a = a z = 5) yield t?£ ~ 0.4 ms and tf£ ~ 0.1 s. 

The additional errors resulting from this finite acqui- 
sition time are easily estimated. We do so here for hard 
spheres (HS). The short time diffusion leads to an error 
Shs = y/D s ((j))t im /2. Here D s (4>) is the volume frac- 
tion (</>) dependent short time diffusion constant for hard 
spheres: D s {4>) = (k B T) / (6nna^HU) with H(d>) < 1 a 
hydrodynamic correction 0, Il36l . 1 1371 . Il38l Il39j |. Note 
that this is an upper bound applicable to HS; for col- 
loids with softer interactions, 5 < 5hs- With a solvent 
viscosity r\ ~ 3 x 10~ 3 Pa-s, and frame rates as above, 
typical values are 82D — 2 nm and 83D — 35 nm for a 
colloid with radius a = 1 /an. For 3D imaging, this error 
exceeds the intrinsic sub-pixel accuracy. Further errors 
due to flow-induced distortion of the particle image are 
estimated by comparing the imaging time t lm with the 
time tf required for the flow to displace the particle over 
its own diameter. For a flow velocity V (in pixels per 
second), tf = 2a/V. We consider the particle signifi- 
cantly distorted if t- im /t{ > 0.1. Hence, using Eqs. 1 1121 
the maximum velocities are: 

= 0.1n/ scan , V™f* = 0.1/a/o,. (3) 

Using the above parameters, a typical maximum velocity 
in 2D is V max ~ 500 ^ms _1 , while for 3D images with 
a/a z = 1, we obtain V rnax ~ 2 /jms" 1 . In both cases, 
further improvement could be achieved by removing the 
distortion prior to locating the parti cles by correlating 
scanned images or lines at different z jlll| . 

C. Tracking algorithms 

Merging particle coordinates from subsequent frames 
into single-particle trajectories is the optimal route to 
analyse co lloid al dynamics. However, in some cases this 
is difficult 140], but quantitative information may still be 
obtained via a lternative methods. For example, Breed- 
veld et al. [H, ll4Ct Il4ll ] measured the distribution of all 
possible displacement vectors in images of tracer particles 
in dense non-Brownian suspensions. From these distri- 
butions, the contribution due to cross correlations, corre- 
sponding to vectors connecting the position of particle i 
in one frame to that of particle j in the next frame, could 
be subtracted. From the resulting autocorrelation part 
of the distribution, the full shear-induced self-diffusion 
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tensor could be obtained. While this method also has 
potential for analysis of coordinate ensembles obtained 
from (confocal) microscopy during flow, we will focus on 
a complete analysis of the particle dynamics via tracking. 

1. Conventional Tracking 

The most widely used algorithm to track particles from 
ensembles of coordinates in consecutive frames is that of 
Crocker and Grier (CG, [l33|). Jt 

is based on the dynam- 
ics of dilute non- interacting colloids. Given the position 
of a particle in a frame and all possible new positions in 
the next frame, within a 'tracking range' Rt of the old 
position, the algorithm chooses the identification with 
the minimum mean squared frame to frame displacement 
(MSFD). The algorithm has been used to analyse parti- 
cle dynamics in a wide variety of 2D and 3D images of 
quiescent systems, see e.g. [TH, Its main limitation is 
that, when the particle motion between frames is exces- 
sive, misidentifications can occur. Such motion can be 
due to diffusion, flow-advection or both. 

In quiescent systems, the ability to track particles be- 
tween frames is limited by experimental constraints such 
as the acquisition rate of the (confocal) microscope, im- 
age dimensionality, particle size and concentration, and 
solvent viscosity. After locating the particles, the rele- 
vant quantity which relates to the tracking performance 
of the CG algorithm is the root mean squared frame- 
to-frame displacement relative to the mean interparti- 
cle distance t. If particles move on average a substan- 
tial fraction of £, then the algorithm starts to misiden- 
tify them. This has recently been quantified [961 ] by 
testing the CG algorithm on computer-generated hard- 
sphere or hard-disk ensembles at different concentrations. 
These tests showed that the algorithm can handle larger 
mean squared displacements (MSDs) for more concen- 
trated systems, since the higher concentration prevents 
particles from coming into close proximity of each other 
between frames. 

The tests were also performed on computer-generated 
data in which additional uniform or non-uniform motion 
was added, to study how far the CG algorithm could 
be pushed beyond its original design parameters. For 
uniform motion, CG tracking was as successful as in the 
quiescent case for small drifts but failed for drifts of the 
order of half the particle-particle separation. For non- 
uniform (linear shear) flows with small strains between 
frames the identification worked correctly, but large non- 
uniform displacements caused major tracking errors. 



2. Correlated Image Tracking 

Some important limitations to particle tracking in 2D 
or 3D images with large drift or non-uniform motion (see 
also [l42j ]) can be overcome by the method of Correlated 
Image Tracking (CIT), described in detail in [96]. The 



main extension compared to conventional (CG) track- 
ing is that, prior to tracking, the time- and position- 
dependent advective motion is obtained from a PIV-type 
correlation analysis as in Section IIII Al This advected 
motion is then subtracted from the raw particle coordi- 
nates, shifting them to a 'locally co-moving' ('CM') ref- 
erence frame where the particles can be tracked with the 
CG algorithm. 

Since particle coordinates are distributed continuously, 
the advection profile AR(x p , y q ,z r ,ti) form the PIV anal- 
ysis is first interpolated to obtain a continuous profile 
AR(i, y, z, ti). Using the latter, the transformation of 
the position r k (ti) = [xk(U),yk(U), 2fe(*i)J of particle k 
in the laboratory frame to its position Ffe(tj) in the CM 
reference frame is: 

i 

r k {U) = r k {U) - J2 AR(r k (t;), t s ), (4) 
j=i 

with AR(rfc(tj), tj) the past motion between frame j and 
j — 1, at the current particle location r k (ti). In the CM 
frame, the average particle motion is essentially zero over 
the entire image. Therefore the classic CG algorithm 
tracks particles successfully in this reference frame. The 
limitations on the CG tracking performance are essen- 
tially the same as discussed in Section IIII C II for qui- 
escent systems. Finally, the resulting trajectories can be 
restored in the laboratory frame by inverting Equation^ 
The main limitation to CIT originates from the fail- 
ure of the image correlation procedure when the relative 
particle motion between frames is excessive, rather than 
from failure due to large absolute shifts between images 
[96| . This apart, CIT has the same limitations of CG 
tracking with significant failure for mean squared frame- 
to-frame displacements in the CM reference of ~ (O.M) 2 
in 2D. Two direct applications of CIT will be reviewed 
in Section ITVAl 

IV. IMAGING OF SYSTEMS UNDER 
DEFORMATION AND FLOW 

Direct imaging remains the most detailed method for 
studying structure and dynamics of particles in flow. 
This section reviews the application of (confocal) imag- 
ing to study flowing concentrated suspensions. We sep- 
arately discuss disordered systems and systems in which 
order is present before or as a result of flow. In each case, 
the focus is on Brownian systems, but we briefly review 
non-Brownian systems at the end of each subsection. 

A. Disordered systems 

1. Slow glassy flows 

A significant motivation for recent developments in 
imaging concentrated colloidal suspensions under flow is 
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to investigate the behaviour of samples that are so dense 
that structural rearrangements are arrested in the quies- 
cent (unsheared) state, i.e. colloidal glasses ||. Eluci- 
dating the mechanisms for the deformation and flow of 
colloidal glasses is currently one of the 'grand challenges' 
in soft matter science. Much insight to this problem can 
be expected to come from the detailed study of model 
systems. The simplest model systems are concentrated 
hard sphere suspensions, which are glassy at volume frac- 
tions 4> ^ 0.58. (This claim has recently been disputed 
[l43j ] : but however this controversy eventually resolves, 
it remains true that for almost all practical rheological 
purposes, amorphous colloidal hard-sphere suspensions 
behave as soft solids for (j> > 0.58.) 

The first time resolved studies of the flow of HS col- 
loidal glasse s at single particle level have been described 
in [96|, [93, Il04| . Both the 3D particle dynamics and 
global flow were investigated in steady shear under var- 
ious boundary conditions usi ng bo th a planar shear cell 
and 'confocal rheoscope' [9?], Il04l | (geometries with uni- 
form stress in the flow gradient direction z). Standard 
fluorescent PMMA-PHS particles (a = 850 nm, = 0.62) 
in a charge-screened, CHB-decalin mixture were used 
for observations at the particle scale, measurements on 
a global scale employed fluorescent tracers in a non- 
fluorescent host suspension. 

We first focus on dynamics observed for rough bound- 
ary conditions. Globally, these highly concentrated 
HS glasses shear non-u niformly, showing rate-dependent 
shear localization [97|, Il04 Il44j : at the lowest applied 
rates 7 a , the local shear rate 7 near (one of) the walls con- 
siderably exceeds ja and vanishes smoothly into the bulk 
where the system remains solid-like. The characteristic 
length, along z, characterizing this shear rate variation 
starts from ~ 25 particle diameters and increases with 
the applied rate/stress, i.e. the wall 'fluidization' propa- 
gates into the sample and eventually leads to fully linear 
flow profile at the largest rate. 

The microscopic dynamics inside a shear band was 
then studied by focussing on a 30 x 30 x 15 /1m 3 volume 
(~ 3000 particles) the bottom layer of which was ~ 10 
particles away from the wall. A series of undistorted 3D 
snapshots of the entire microstructure was acquired, from 
which 3D particle positions and dynamics were obtained 
as described in Sees. IIII Bl and IIII CI Particle tracking 
using the new CIT method allowed to study the flow for 
local shear rates up to 7 ~ 0.05 s _1 . The results are 
shown in Figure [5J where (a) presents the advected mo- 
tion AX(z, At) (from image correlation) for At = 16 s 
and a shear rate 0.02 s _1 96]. Local velocity profiles 
such as these are linear on this scale and extrapolate to 
zero within the resolution, showing the absence of slip 
in this case. Figure [SJb) shows the distributions P(Ax) 
and P(Ay), which are frame-to- frame particle displace- 
ments obtained from CIT. The laboratory frame data 
for the displacements in the velocity direction (ir) show 
a large contribution from the advected motion. To fo- 
cus on the more interesting non-affine part of the dis- 




FIG. 8: Reprinted from [96j. 3D analysis of a sheared 
glass at 7 = 0.019 s _1 . (a) (•) accumulated displace- 
ment AX a (z, At) (see Figure |6]) from image correlation over 
At — 16 s (4 frames). Line: linear fit giving an accumulated 
strain dAX a /dz = 0.28. (b) Distribution of frame-to-frame 
displacements P(Ax) and P(Ay) from Correlated Image 
Tracking with coordinates restored in the laboratory frame. 
Also shown is P(Ax) of the non-affine ^-displacements, us- 
ing Equation [5] and 7 = 0.019 s -1 . (c) the (non-affine) 
MSD in the three directions. Line: (Ay 2 (At)) = 2Dt with 
D = 5.4 x 10~ 3 /WV 1 . 

placements, AS, the ^-dependent advected motion was 
subtracted via: 

Ax(t) = x(t) - x(0) - 7 / z(t')dt', (5) 
Jo 

with 7 measured form the data. The results for P(Ax), 
Figure [FJb), showed that non-affine motion for x and 
y were very similar. Analysis of the full 3D dynamics 
further revealed that it is indeed nearly isotropic in all 
directions, as shown by the long time behaviour of the 
mean squared displacements for x, y and z in FigurelSJc). 
The linearity of these curves at long times also indicates 
that, after a sufficient number of shear induced cage- 
breaking events, particle dynamics in the system becomes 
diffusive, with long-time diffusion coefficients varying by 
< 20% in the three directions. 

These experiments give information on individual plas- 
tic cage-breaking events as well as rate-dependent struc- 
tural relaxation. The first phenomenon is illustrated by 
the non-affine motion of a single particle in the x — y 
plane, Figure [9^a) . Intermittent jumps, reflecting 'plas- 
tic' breaking of the particle cages, are observed between 
periods of 'rattling' in which the cage is deformed. An 
average accumulated strain of ~ 10% is found between 
these events, in reasonable agreement with the yield 
str ain o btained from bulk rheology or light scattering 
[1, I XOOj ] . The changes in the particles' local environ- 
ment during shear were also studied via the chan ge in 
the bond order parameters Q l 6 m for each particle i jl45| 
over a time dt, measured by the quantity C\{t, dt), where 
Cq = reflects no change. Figure Efc)-(e) show snap- 
shots of the sheared microstructure at two times, where 
the grey-scale is proportional to C|(t, dt). In Figure^c), 
the changes over 6% accumulated strain are shown. Clus- 
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FIG. 10: Reprinted from [l 52t [ . with permission from AAAS. 
(A) Shear-induced displacements of particles in a colloidal 
packing between z = (the cover slide) and z = 23 /im after 
50 minutes of shear (3% average accumulated strain). (B) 
Strain distribution and shear transformation zones in a 7 /im 
thick section in the shear-displacement gradient plane. Parti- 
cle color indicates the value of the local shear strain e yz (see 
color scale) , accumulated over 50 minutes of shear. The arrow 
indicates a shear transformation zone. 



FIG. 9: Partly reprinted from [1041]. (a) Trajectory of a sin- 
gle particle over 800 s in the velocity-vorticity plane in the 
co-moving frame as defined in Equation [5] for a sheared glass 
at a local shear rate 7 = 9.3 • 10~ 4 s _1 . (b) Selected inco- 
herent scattering functions F s (Q m ,t), with 7 increasing from 
right to left. The dashed line schematises initial relaxation. 
Inset: data collapse using fs(Qm,t/T a ), with the line show- 
ing that f s tx exp(— t/r a ). (c-f) Snapshots of the structure 
for 7 = 0.0015s" 1 , (c) at to = 80s after start of imaging. The 
grey-scale measures the change in local environment of each 
particle (oc Ce(to, dt)) over the past dt = 40s. (d) at to = 120s 
and dt — 80s. (e) at to = 120s and dt = 40s. Red lines show 
local deformations, arrows mark the shear direction, (f) A 
quiescent glass, with dt = 200s. 



ters of strong rearrangements, termed 'Shear Transfor- 
mation Zones' (STZs) in earlier theoretical studies jl46| . 
can be observed, where the local strain appears much 
higher than elsewhere. Over the following 6% accumu- 
lated strain, Figure [9](e) , STZs appear in different loca- 
tions while the earlier ones remain essentially locked. 

The average structural relaxation was examined 
via the incoherent scattering function, F s (Q,t) = 
(cos(Q[yi(to+t)-yi(to)])} iito (using Q = Q m ~ 3.8a- 1 ). 
As shown in Figure [9jb), the long-time dynamics is es- 
sentially frozen at rest (7 = 0). At short times and small 
shear rates, F s exhibits a plateau corresponding to caging 
at small accumulated strain. At longer times, F s decays 
to zero due to repeated cage-breaking events in line with 
the diffusive dynamics in Figure (He). This decay accel- 
erates strongly on increasing rate. When scaling time 
by the characteristic (a) relaxation time r Q , F s collapsed 
onto a single exponential curve (see the inset), which di- 
rectly confirmed the theoretically predicted time-shear 
superposition principle. 

The rate dependence of the inverse relaxation time and 
diffusivity D exhibited a non-trivial scaling: D ~ l/r a ~ 
7 8 . Such 'power-law fluid' scaling contrasts the yield- 
stress behaviour observed in the glob al rh eology or pre- 
dictions from Mode Coupling Theory jl47| | . (For further 



discussion of these data in the context of Mode Coupling 
Theory, see the article by Fuchs in this issue.) However, 
these observations are consistent w ith a sto chastic non- 
linear Langevin equation treatment jl48lll4t^ |. In this ap- 
proach to hard-sphere dynamics, particles become more 
deeply trapped in effective free-energy (or, equivalently, 
entropic) wells as the concentration increases beyond a 
certain critical value; the effect of applied stress is to 
lower the entropic barri er that particles have to surmount 
in order to escape |150| : the predicted shear-induced re- 
laxation time indeed shows power-law scaling with shear 
rate s in an intermediate regime with exponents close to 



0.8 [151 



Subsequently, Schall et al. studied the behaviour of 
dense collo idal packings under very small shear deforma- 
tions [l52]. The packing was formed by sedimentation of 
silica-spheres (a = 0.75 /im) in a lower density solvent, 
so that the volume fraction varied as function of height. 
The shear geometry consisted of a movable, coated bot- 
tom cover slide, while a metal grid positioned on top 
confined the packing within a layer of 42 /zm (~ 30 parti- 
cle diameters) . The analysis focussed of the microscopic 
strain variations resulting from either thermal fluctua- 
tions or applied strain. From the thermal distribution of 
strain energies, an estimate for the shear modulus was 
obtained, but volume fraction gradients complicate the 
analysis. Large accumulated local strains were observed 
at long times in the quiescent sample, indicating aging. 
The behaviour under application of shear is shown in Fig- 
ure [TO] In this work, the total accumulated strain was 
~ 3% on average, Figure [TWa). well below the t ypica l 
bulk yield strain in colloidal hard sphere glasses [8|, llOOj ] . 
Thus, it is probable that the data represent creep rather 
than yielding and flow. 

Figure fTW b) shows a spatial map of the shear strain e yz 
in a slice normal to the shear gradient. The authors iden- 
tified the localized regions with high strain with STZs, 
the first published images of STZs in colloidal glasses. 
Interestingly, the highlighted region revealed a four-fold 
symmetric strain field just after its (earlier) formation, 
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FIG. 11: Reprinted from [97]]- (a) Measured stress cr m versus 
applied rate -y a for particles with a = 138 nm, at different 
volume fractions for coated (■,*,▲) and un-coated geome- 
tries (DjOjA). Dashed line: linearly viscous behaviour. The 
full line is a flow curve predicted using the model in [97j ] . 
Dotted line: applied rate at which yielding starts at the cone- 
edge, (b) Density profile P{z) (□), from 3D-imaging of the 
a = 652 nm system at <f> = 0.63, during slip at r = 2.5 mm 
and % = 0.01 s _1 . (•) Corresponding velocity profile from 
particle tracking, showing full plug flow. Inset: Slice of one 
of the 3D-images, showing the first colloid layer. Scale bar: 
10 /im. (c) Velocity profiles for the (f> — 0.59 data in (a), in 
units of the cone velocity t> CO nc = 7a Or, as function of the 
normalized hight z/z sap — z/Or, for various applied rates 
7a at r = 2.5 mm. The arrow marks the slip velocity for 
7a = 30 s - ; the dashed line is the behaviour without slip 
v = 7a z as observed for <j> < 0g- 



This is confirmed by imaging (Figure ITlTc). dashed line) 
and also agrees with earli er ve locity profile measurements 
for dense colloidal liquids [lll| . However, for the colloidal 
glass, a slip branch developed in the flow curve at small 
rates. This branch was described by a Bingham-form, as- 
sociated with solid-body motion of the suspension along 
the glass plate. This 'solid' structure during slip was 
shown to extend down to the first particle layer, Fig- 
urc lllf b) and the inset, giving microscopic insight in the 
'Bingham' slip-parameters. The transition from pure slip 
at low rates to yielding at high rates was also probed, see 
Figure UlT c). This transition depended on the local gap, 
i.e. the radial position r, illustrating the possibility of 
non-uniform (r-dependent) stress in a cone-plate. 

Interestingly, this slip behaviour of hard-sphere glasses 
is different in nature from that found earlier by Meeker et 
al. in jammed systems of emulsion droplets [60(. There, 
a non-linear clasto-hydrodynamic lubrication model, ap- 
propriate for deformablc particles, could quantitatively 
account for their observations. It therefore appears that 
while slip is ubiquitous, the mechanism for its occurrence 
can be highly system dependent. 



consistent with theoretical predictions for a single plas- 
tic event. Schall et al. further identified the STZ-cores 
and, using somewhat uncontrolled assumptions, calcu- 
lated the STZ formation energy, activation volume and 
the activation energy. From this it was concluded that 
STZ formation was mainly thermally activated, although 
some strain assistance was involved. 

Images such as those shown in Figure fTOlb) and[9jc)- 
(e) show that direct imaging is a very powerful method 
for studying inhomogeneous local responses to applied 
stress. In particular, the discovery of STZs in model 
colloids means that the physics of these systems may have 
relevance beyond soft matter physics, e.g. i n the study 
of large-scale deformation of metallic glasses [1531 ] . 

Particulate suspensions and other concentrated soft 
materials (e.g. emulsions) often exhibit wall slip dur- 
ing flow along smooth boundaries. This is important in 
practical applications, and strongly relevant when inter- 
preting bulk flow measurements. Slip in particulate sus- 
pensions has been seen mostly for high Peclet num bers 
(non-Brownian suspensions), both in solid- |l54ll55( and 
liquid- like [l56j | systems. In a recent study, Ballesta et 
al. addressed the effect of Brownian motion and the glass 
transition on slip in hard-sphere colloids using confocal 
imaging and simultaneous rhcology (9?| (Section III C 2[) . 
The main results are shown in Figure 1111 The rheol- 
ogy for rough walls in Figure [TlT a) shows the traditional 
change from a shear thinning fluid (4> = 52%) to a yield- 
ing solid (</> = 5 9%) associated with a glass transition 
at (f> ~ 58% 0, Il57| |. For a smooth wall and index- 
matching solvent (which prevents van der Waals attrac- 
tion and sticking of particles to the glass) the results for 
a dense liquid remained unchanged, indicating no slip. 



2. Fast channel flow of fluids and pastes 

A different regime of flow in hard-sphere suspensions, 
cither in dense fluids or glasses, occurs when the shear 
rate becomes comparable to the rate for thermal relax- 
ation of the particles within their cages, which is of the 
order of the inverse Brownian time. In these situations 
jamming or shear thickening may start to occur, while in 
geometries with non-uniform stress, pronounced shear- 
induced migration may take place. Here we describe re- 
cent studies of imaging in this regime, focussing on the 
specific geometry of capillaries and micro-channels. 

The effects of shear-induced migration in intermediate 
volume fraction hard-s pher e suspensions {<f> < 0.35 ) wa s 
studied by Frank et al. [1161 ] and Semwogerere et al. [11 71 ] . 
In these experiments, confocal imaging was used to mea- 
sure both velocity and concentration profiles, <f>(y), of 
colloids across 50/xmx500/im glass channels. Flow ve- 
locities up to 8 mms -1 were probed, corresponding to 
Peclet numbers up to ~ 5000. The gradient in shear 
rate j(y) causes particles to migrate from the bound- 
aries (with large 7) to the center of the capillary (small 
7), resulting in development of a steady state concen- 
tration profile beyond a certain 'entrance length'. Frank 
et al. found profiles (transverse to the flow) where the 
concentration in the center increased considerably with 
the average volume fraction or on increasing flow rate 
(Peclet number, Pc). The data were interpreted via a 
model which included a shear rate and local volume frac- 
tion dependence of normal stresses in the sample. The 
rapid rise of these stresses with <fi and Pe was responsi- 
ble for the observed behaviour. Semwogerere et al. |l 17j ] 
studied the entrance length over which the flow became 
fully developed and found that it increases strongly with 
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FIG. 12: (a) Velocity profile of a 63.5% volume fraction sus- 
pension of PMMA hard-spheres (radius a = 1.3/im) in a 
square, glass micro-channel (side 2b — 50/im). The velocity 
(V) is normalized by the centerline velocity (V c ). The profile 
is overlaid onto a confocal image of the suspension in the chan- 
nel. The velocity profile consists of a central, unsheared plug 
and of lateral shear zones, whose width is highlighted by the 
horizontal blue lines. The right side of the image shows two 
examples of particles trajectories in the plug and in the shear 
zone respectively. The x and y scales of the trajectories are 
differ ent a nd the colloids are not drawn in scale, (b) Reprinted 
from [129([ . Normalised velocity profiles at z = 17 /im versus 
y/b e g (b e fi = b — a) for a wide range of flow rate. For clarity, 
fits to the data rather than the raw data are shown; the error 
bars show the spread in the measurements. Inset: normalized 
transverse fluctuations of individual particles in the shear zone 
for the two boundary conditions, versus normalized time t/rs 
or local accumulated strain t ■ d{V(y)) jdy = t ■ 7(1/), top axis. 



Pe, in marked contrast to n on-B rownian flows for which 
it is flow-rate independent [158]. The entrance length 
also decreased with increasing volume fraction and the 
data could be successfully described within their model. 

Few imaging studies have bee n car ried out in higher 
concentration suspensions. Haw jl!2| | observed intrigu- 
ing behaviour in traditional microscopy studies of PMMA 
suspensions (4> > 55%) sucked into mm-sized capillaries 
with a syringe. Strikingly, the volume fraction of the sus- 
pension in the capillary was as much as 5% lower than 
the initial concentration, depending on the particle size. 
Such unexpected "self- filtration" was attributed to jam- 
ming of the particles at the capillary entrance, leading 
to higher flow rate of the pure solvent under the applied 
pressure difference. 

Motivated by these experiments, Isa et al. conducted 
further studies of flows in similar geometries at th e sin- 
gle particle level using confocal microscopy [129]. The 
system consisted of a hard-sphere suspension (PMMA 



spheres, radius 1.3 ±0.1 fim) at nearly random close pack- 
ing, 'a paste', in a 20-particle-wide square capillary. The 
motion of individual colloids was tracked via Correlated 
Image Tracking and velocity profiles were measured in 
channels with both smooth and rough walls. Despite the 
colloidal nature of t he suspen sion, significant similarities 
with granular flow [l59l Il60l | were found. 

The bulk flow curve of the system studied by Isa et 
al. fitted a Herschel-Bulkley (HB) form for a yield-stress 
fluid [1] at small to moderate flow rates. The velocity 
profile predicted from a HB constitutive relation consists 
of a central, unsheared "plug" and shear zones adjacent 
to the channel walls, which fits qualitatively with what 
was observed, Figure [T2Ta) . At the same time, however, 
the size of the central plug is expected to decrease with in- 
creasing flow rate, vanishing beyond some critical veloc- 
ity above which the yield-str ess fluid is shear-melted ev- 
erywhere (see e.g. [l05l ll6l| V But Isa et al. found that 
width of the central plug was independent of applied flow 
rate, Figure [T27b). This behaviour is analogo us to tha t 
observed in the pipe flow of dense dry grains jl6dl Il62l ] . 
The data for the colloidal flow could be captured by a 
theory of stress fluctuations origi nally dev eloped for the 
chute flow of dry granular media jl62L Il63j . Presumably, 
this model can be successfully applied to nearly-close- 
packed colloidal systems because at high flow rates, inter- 
particle contacts rather than Brownian motion dominate. 
Direct imaging showed indeed that the trajectory of a 
particle in the shear zone was determined by 'collisions' 
with neighbouring particles, Figure [T2j Moreover, the 
presence of rough boundaries enhanced such fluctuations, 
causing wider shear zones. 

In more recent work on the same system [ll3j . us- 
ing smooth walls, it was shown that confinement can in- 
duce instabilities in the flow. Upon decreasing the width 
of the channels from ~ 40 to ~ 20 particle diameters, 
the flow developed oscillations above a threshold applied 
flow rate. Such oscillations consisted of cyclical jamming 
and un-jamming of the suspended particles and l ed t o 
filtration effects similar to those reported by Haw [l!2| . 
Single-particle imaging was used to demonstrate the pres- 
ence of a concentration profile across the channel, which 
was well correlated with the local velocity profile. 



3. Non-Brownian suspensions 

Other recent imaging experiments address the flow of 
non-Brownian suspensions (i.e. Pe — > 00) at the singl e 
particle level. In a series of papers [13, IH, Il40l Il4l| . 
Breedveld et al. investigated the steady shear-induced 
self diffusion of non-Brownian spheres. They used 
PMMA particles (radius a ~ 45 /itm) in a RI and density- 
matching mixture (water, zinc chloride, and Triton X- 
100, viscosity ~ 3.4 Pa s) and studied suspensions with 
volume fractions ranging from 20% to 50%. By studying 
the correlated motion of a small fraction of the particles 
which had been dyed (Section IIII C|) . the self and 'off- 
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diagonal' mean squared displacements (MSDs) could be 
extracted without the use of explicit particle tracking. 

For these non-Brownian suspensions, the long-time 
shear-induced MSD depend on accumulated strain jAt 
only. The 'long-time' diffusion for ^yAt > 1 is due to the 
chaotic nature of multiple particle collisions. The asso- 
ciated self-diffusion constants D xx , D yy , D zz (in the ve- 
locity, vorticity and gradient direction respectively) were 
found to be anisotropic with D zz /D yy < 2, while D xx 
was almost an order of magnitude larger. The data 
were in qualitative agreement with hydrodynamic theo- 
ries and simulations, although quantitative discrepancies 
remained. In c ontra st, for the slowly sheared glasses in 
Section ELO p^3| . which exhibited a non-trivial rate 
dependence, D yy /(a 2 j) scales as 7~ ' 2 , ranging from 0.7 
at low rates to ~ 0.3 at the largest rates, exceeding 
the result for non-Brownian suspension. A l ast n otable 
observation in the non-Brownian system in [l4dj | was a 
regime of reduced diffusion for small accumulated strains 
'j At ~ Ajcage <C 1. While its nature remains unclear, it 
was shown that the typical strain A7 CQge roughly corre- 
sponded to the affine deformation S/a at which particles 
come in direct contact (on average), which varies with 
volume fraction as (c/yep/^) 1 / 3 ~ 1; an d vanishes as the 
average surface-surface separation between particles on 
approaching close packing. 

In a very recent study by Wang et al. [l09t | , the three- 
dimensional motion and diffusion of heavy tracer par- 
ticles in an immersed, density matched granular pack- 
ing {(j) ~ 0.6) in a cylindrical Couette-cell was stu died . 
Among many other remarkable results, the data in [l09j ] 
showed direct evidence that D ex 7 in non-Brownian sys- 
tems. Interestingly, the value for shear diffusion in the 
vorticity direction, Dj (7a 2 ) ~ 0.1 was close to the values 
found by Breedveld for slightly lower density. 



Another study on non-Brownian suspensions 164| fo- 
cussed on a transition from reversible to irreversible par- 
ticle dynamics in oscillatory shear. While the creeping 
(Stokes) flow equations for simple shear flow are in princi- 
ple reversible, particle roughness, three particle collisions 
or repulsive forces (ignoring Brownian motion) may ren- 
der the dynamics irreversible. Pine et al. showed that 
there exists a volume-fraction-dependent critical strain 
amplitude beyond wh ich t he dynamics become diffusive. 
In a follow-up study [l65j . they demonstrated that the 
state below the critical strain is in fact an 'absorbing 
state', where the particles self-organize into a structure 
where collisions are avoided. They also demonstrated 
an increase in the dynamic viscosity of the suspension 
on crossing the critical strain amplitu de, a s was also ob- 
served previously by Breedveld et al. jl40l |. An interest- 
ing question is whether any relation exists between the 
absorbing states measured in these non-Brownian sys- 
tems and the behaviour below the yield strain in colloidal 
glasses. 

To date, flow studies of Brownian and non-Brownian 
systems have largely been conducted independently. It 
is clear that much can be gained by detailed comparison 



of such studies [166] and that single-par ticle imag ing can 
provide crucial insights in this exercise [l29l Il67j . 



B. Ordered systems and ordering under flow 

Hard spheres colloids undergo an entropy driven fluid- 
crystal phase transition when the suspension's volume 
fraction reaches <j)f = 0.494. This behaviour was 
first confirmed experimentally using quasi-monodisperse 
hard-s phere - like PMMA colloids [6J. Static light scat- 
tering [l68j | shows (and later, conventional microscopy 
[sTj ] confirms) that the resulting colloidal crystals con- 
sist of hexagonally-packed layers more or less randomly 
stacked on top of each other; i.e., these crystals have a 
mixture of face-centered cubic (fee) and hexagonal close 
pack ed (hep) structures. Starting with a paper in 1998 
[l69j |. van Megen and co-workers have published a se- 
ries of studies in the kinetics of crystallization in PMMA 
suspensions using time-resolved dynamic and static light 
scattering, inv estigating in particular the effect of poly- 
dispersity [lTOj ] in some detail. 

Gasser and co-workers [H| pioneered the use of con- 
focal microscopy to study the kinetics of colloidal crys- 
tallization. Using buoyancy-matched PMMA particles, 
they observed nucleation and growth of colloidal crys- 
tals with single-particle resolution. These authors iden- 
tified the size of critical nuclei, 60-100 part icles, in rough 
agreement with computer simulations [17l| , and to mea- 
sure nucleation rates and the average surface tension of 
the crystal-liquid interface. For completeness we men- 
tion that Palberg and coworkers have pioneered the use 
of non-confocal microscopy metho ds to image ordering in 
highly-deionized suspensions [13, Il72| , which crystallise 
at very low volume fractions (10~ 3 or lower). 

The relative mechanical weakness of colloidal crystal 
compared to their atomic and molecular counterparts 
means that modest shear will have large, non-linear ef- 
fects on crystallisation as well as crystal structure. In 
the case of hard-sphere-like PMMA colloid s, th ese effects 
have been studied using light scattering jl73| since al- 
most immediately after the discovery of crystallisation in 
this system Q . Much more recently, confocal microscopy 
has been used to image shear effect in PMMA colloids. 
Below, we focus on reviewing such studies, but we note 
that imaging studies of crystallisation in other colloids 
have also been performed (see, e.g., (13)- To date, there 
have been few imaging studies of shear effects on order- 
ing and ordered states in non-Brownian suspensions. We 
review briefly one series of such work at the end of this 
section. 



1. Hard-sphere colloids 

It is known from static li ght scattering [l73l |. diffus- 
ing wave echo spectroscopy [86] and conventional opti- 
cal microscopy (85j that the application of oscillatory 
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FIG. 13: Reprinted with permission from [9^]. Copyright 
[2006], American Institute of Physics. Confocal images (top 
row, with insets showing the Fourier transform) and recon- 
structed 3D stacks (bottom row) showing the time dependent 
effect of oscillatory shear (frequency = 3 Hz, strain amplitude 
= 1) on a 48% volume fraction PMMA suspension. Time pro- 
gresses from left to right: quiescent (24 hours after prepara- 
tion, A), after 2 minutes (B), and after 5 minutes (C) of shear. 
In (B) and (C), unit vectors (e) in the the velocity (v), vor- 
ticity (id) and shear gradient (V) directions are shown. The 
light grey particles in the 3D renderings (bottom row) rep- 
resent those with a high degree of local crystalline ordering. 



shear to dense colloidal fluids or glas ses can drive crys- 



tallisation. Rheological studies [174j show that shear- 
crystallised samples have lower elastic and viscous mod- 
uli than their glassy counterparts at the same volume 
fraction, in agreement with the enhanced effective 'cage' 
volume in the ordered case, leading to smaller entropic 
stiffness. 

A detailed confocal imaging study of shear-induced 
crystallisation in hard-sphere colloids was performed re- 
cently by Solomon and Solomon [94j . They imaged crys- 
tals formed under an oscillatory shear field at a particle 
volume fraction of 52% in slightly-charged PMMA par- 
ticles (diameter 1.15 /im, polydispersity 4%) in a planar 
shear cell (Section III C 1|) . Particles were identified using 
the conventional Crocker and Grier algorithm but, due 
to the low acquisition rate (0.6 frames/second), imaging 
was only possible after the cess ation of shear. Consis- 
tent with previous work [85l Il73| | , they observed ordering 
of the initially amorphous suspension into close-packed 
planes parallel to the shearing surface, Figure [T51 Upon 
increasing the amplitude of the applied oscillatory strain, 
7, the close-packed direction of these planes was observed 
to shift from an orientation parallel to the vorticity direc- 
tion to one parallel to the flow direction, and the qual- 
ity of the layer ordering decreased. In addition, they 
studied shear-induced stacking faults and reported their 
three dimensional structure. For large STRAIN ampli- 
tudes (7 = 3), ordering in the flow- vorticity plane only 
persisted for 5 to 10 layers, while in the gradient direc- 
tion the crystal consists of alternating sequences of fee 
and hep layers. 

Once formed, whether it be under quiescent or sheared 
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FIG. 14: Reprinted with permission from [llll ]. copyright 
[2004], Institute of Physics Publishing (a) Superposition of 
particle tracks in the zero-velocity plane plus one adjacent 
plane in a sheared colloidal crystal, showing zig-zag motion. 
The image size is 18.75 fim x 18.75 fj,m (256 x 256 pixels). 
The colloids have a diameter of 1.50 /^m. (b) Measured shear 
rate versus applied shear rate for the colloidal crystal (points) , 
the continuous line shows equality of these two rates. 



conditions, colloidal crystals show a complex response to 
applied stress. In particular, they exhibit banding under 
simple shear. We have already discussed flow localiza- 
tion in colloidal glasses in Section TlV A II The first evi- 
dence for the occurrence of this phenomenon in ordered 
structures was reported in a study of hard-spher e col - 
loidal crystals under steady shear by Derks et al. 
They studied the velocity profiles and single-particle dy- 
namics via correlation techniques and particle tracking 
in a cone-plate system (Section III C 2[) around the zero- 
velocity plane. Since particles were locked in their lat- 
tice positions and their motion was limited, they could 
be tracked using the conventional Crocker-Grier algo- 
rithm. The close-packed direction in hexagonally-ordered 
planes was observed to be parallel to the velocity, in 
agreement with high strain amplitude oscillatory shear 
results [H, ■ The study also showed direct evidence 
that colloids in adjacent layers move in a zig-zag fashion, 
Figure [TJJa), as previously inferred from light scatter- 
ing jl75l ] and conventional imaging [85| . The velocity 
profile measured over the accessible z range were linear, 
but the measured shear rates were much larger than the 
applied rates, Figure [F4T b). indicating considerable shea r 
banding in the sample. In a more recent study [l76j . 
Derks et al. analysed in more detail the particle dynam- 
ics in the flow-vorticity plane of a sheared crystal. They 
evidenced increasing fluctuations in the crystal layers on 
increasing rate, eventually leading to shear melting when 
the short-time mean square displacements due to shear 
attained a 'Lindemann' value of ~ 12% of the particle 
separation. 

More direct evidence of shear banding in colloidal crys- 
tals und er os cillatory shear was found by Cohen and co- 
workers [l02t | . In this confocal microscopy study, the os- 
cillatory motion across the gap was measured at various 
applied amplitudes and frequencies. For low applied de- 
formation, 7a PP , the authors claim that the crystal is lin- 
early strained at all frequencies, /. At higher 7app, the 
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crystal yields, and two different regimes are observed. At 
low frequencies, a high-shear band appears close to the 
upper (static) plate. As the applied strain and frequency 
increase, the band becomes larger until the whole gap is 
sheared at / > 3 Hz. Moreover, large slip is observed 
in all cases, with the measured strain always below the 
applied one. Cohen et al. proposed a model based on the 
presence of two coexisting, linearly-responding phases to 
explain their observations. This model is attractive for 
its simplicity, but questions remain, in particular con- 
cerning the lack of bulk strain for small-amplitude, low- 
frequency data and some significant deviations between 
the modelled and measured profiles. 

Simple shear is, of course, a highly idealised deforma- 
tion geometry. Motivated partly by the possibility of 
shedding light on the deformation, fatigue and fracture 
of metals in more 'real-life' situations, Schall and co- 
workers carried out a confocal imaging study of a more 
complex form of deformation: indentation in defect-free 
colloidal crystals grown by sl owly sedimenting silica par- 
ticles on templated surfaces |l77| . Indentation (or nano- 
indentation in the case of metallic crystals) consists in 
pressing a sharp tip (diamond for metals, and a needle 
for colloidal crystals) against a surface to probe its me- 
chanical properties by relating the applied pressure and 
the measured penetration. It is also a method for the 
controlled introduction of defects into a crystal. 

Schall et al. indented colloidal crystals using a needle 
with an almost hemispherical tip of diameter 40/im, in- 
ducing a strain field in which the maximum shear strain 
lies well below the contact surface. The tip diameter, 
particle radius and crystal thickness in their experiments 
were chosen to be similar to parameters in typical metal- 
lic nano-indentation experiments. The authors discussed 
their observations using a model that addresses the role 
played by thermal fluctuations in the nucleation and 
growth of dislocations. 

We have already discussed confinement effects in the 
channel flow of colloidal glasses. Such effects are also seen 
in hard-sphere colloidal crys tals sheared between parallel 
plates. Cohen et al. 103] found that when the plate 
separation was smaller than 11 particle diameters, com- 
mensurability effects became dominant, with the emer- 
gence of new crystalline orderings. In particular, the col- 
loids organise into "z-buckled" layers which show up in 
xy slices as one, two or three particle strips separated by 
fluid bands, Figure 1151 By comparing osmotic pressure 
and viscous stresses in the different particle configura- 
tions, the cross-over from buckled to non-buckled states 
could be accurately predicted. 




FIG. 15: Reprinted with permission from 103], copyright 
2004 by the American Physical Society. Structure of a 
sheared suspension (diameter 1.42 /im) with strain ampli- 
tude 70 = 0.38, frequency f — SO Hz and 4> = 0.61. The 
plate moves in the x direction, (a) Confocal micrograph of a 
sheared suspension forming hep layers when the gap between 
the plates D — 80 /im. (b)-(e) images of the suspension in the 
buckled state. The gap is set slightly below the height com- 
mensurate with confinement of four flat hep layers, (b) an 
xy image slice of the suspension near the upper plate, (c)-(e) 
show slices that are, respectively, 1.3, 2.6, and 3.9 fim below 
the slice in (b). 
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FIG. 16: Reprinted with permission from 108], copyright 
2003 by the American Physical Society, (a) Time evolution of 
the stru cture of an index matched suspension of glass beads 
in |l08j ]. imaged in the flow-gradient plane after the start of 
shearing at t — 0. The observed layering is parallel to the flow 
direction, (b) Fractional volume compaction vs time, based 
on the change of height of the upper driving boundary. The 
actual change in microns is indicated at the right. 



sheet imaging was employed to study the dynamics of 
confined glass sphere packings (diameter 600 /im) in an 
index-matching (non-density matching) solvent in an an- 
nular shear cell with load applied from the top shearing 
surface. Using simultaneous measurements of the volume 
(and torque in [l79j |). they discovered a transition from 
a disordered packing to an ordered packing under steady 
shear, Figure [TuT a). accompanied by compaction steps, 
Figure [TBT b). and reduction of the shear force. They 
observed strongly non-linear shear profiles with a rapid 
reduction of the rate away from the top driving plate, 
varying significantly with the order of the packing. Due 
to these non-linear profiles, the delay observed before the 
ordering transition was strongly dependent on the layer 
thickness for fixed boundary speed, since global order- 
ing requires an accumulated strain A7 ~ 1 through the 
entire sample. 



2. Non-Brownian suspensions 

Shear induced structuring has bee n pre dicted in the 
flow of non-Brownian suspension (e.g. [178j ]). but to date 
there have been few imaging studi es. In a series of ex- 
periments by Tsai and co-workers jl08l . Il79l [l80], laser 



V. CONCLUSION AND OUTLOOK 

Recent advances in imaging and data analysis tech- 
niques have enabled time resolved imaging and tracking 
of individual particles in colloidal systems under flow. 
We have reviewed these techniques, as well as examples 
of their application, primarily on concentrated systems 
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of hard-sphere colloids. Despite what has been achieved, 
we believe that this area is only in its infancy. Results 
already obtained raise many more questions than they 
answer. Many phenomena and systems remain little ex- 
plored, and existing results invite further investigation 
using imaging techniques. For example, computer simu- 
lations of crystal nucleation [l8ll | as well as proposed ex- 
planations of complex yielding behaviour in glassy states 
of systems of particles with short-range attraction [iolflU 
can both be directly tested using such methods. In other 
cases, the use of imaging may throw light on well-known 
and long-studied phenomena that remain incompletely 
understood, e.g. shear thickening of concentrated col- 
loidal suspensions The detailed comparative study 
of distinct but related systems using imaging techniques, 
e.g. Brownian and non-Brownian suspensions, or partic- 
ulate suspensions and emulsions, should also be a fruitful 
area for exploration. 

Further technical developments will enhance the power 
of the methodology. Faster scanning methods will ob- 
viously extend the upper limit of flow rates that can 
be st udied . The development of algorithms for identi- 
fying 182] and tracking particles in polydispcr.se systems 
will significantly extend the range of systems that can be 
studied - there has been no quantitative study to date 
of polydisperse suspensions under flow by imaging. Fur- 
ther use of confocal imaging combined with (simultane- 
ous) bulk rheometry will give a wealth of information 
on the relation between micro-scale dynamics and bulk 
flow and provide a crucial test ground for various theo- 
ries. The coupling of three-dimensional imaging with mi- 



crorheological techniques utilising optical trapping ('laser 
tweezers' [183J) should also yield many new discoveries 
and insights; confocal microscopy in two dimensions has 
already been used to explore the structure of the 'wake' 
downstream from a probe p articl e being dragged through 
a concentrated suspension 184]. Studying flow in mi- 
crofluidic geometries using time-resolved single-particle 
imaging is of direct relevance to many emergent lab-on- 
a-chip type applications. And the parallel use of optical 
microscopy and other imaging techniques (perhaps espe- 
cially NMR) to study the same flow system should open 
up many possibilities. 

Finally, the example applications reviewed here pertain 
mostly to well-characterised, model particles. Confocal 
microscopy can fruitfully be used to image the (static) 
micr ostructure of many 'real-life' systems, such as foods 
[l85j |. To date, there have been few con focal microscopy 
and rheological studies of such systems |llOj | . Given the 
importance of rheological properties in the processing 
and utilisation of a wide range of industrial products 
(controlling the 'mouth feci' of food is but one example), 
we may expect increasing application of this methodol- 
ogy to applied colloidal systems. 
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